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A method for the experimental reconstruction of the quan- 
tum state of motion for a single trapped ion is proposed. It is 
based on the measurement of the ground state population of 
the trap after a sudden change of the trapping potential. In 
particular, we show how the Q{a) function and the quadra- 
ture distribution P{x, 6) can be measured directly. In an ex- 
ample we demonstrate the principle and analyze the sensibil- 
ity of the reconstruction process to experimental uncertain- 
ties as well as to finite grid limitations. Our method is not 
restricted to the Lamb-Dicke Limit and works in one or more 
dimensions. 

PACS Nos. 42.50.-P, 42.50Vk 

The central entity of quantum physics is the density 
operator p. It contains all measurable information about 
the state of a system that can be obtained according to 
the principles of quantum physics. Recent theoretical 
advances established constructive procedures to recover 
the full information about the state of a system from the 
repeated measurement of a complete set of observables. 
From the experimentally detected probabilities 

Q{a) = -Tr[p\a){a\] = -{a\p\a), (la) 
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P{x, 9) ^Tiip\x,e){x,9\]^ {x, 9\p\x,9), (lb) 

one can determine the state p uniquely. Here \a) de- 
notes a coherent state and \x, 9) is a quadrature eigen- 
state. So far, the underlying theory has been developed 
for finite dimensional discrete systems [|l|,H, like spin- or 
angular momentum states, as well as for continuous sys- 
tems. This approach is generally referred to as phase 
space tomography 

One of the most important and beautiful applications 
has been the tomographic measurement of the Wigner 
function for a single mode of the electro-magnetic radia- 
tion field by Raymer et.al. following the proposal by 
Vogel and Risken ||^ . It is based on a measurement of the 
quadrature distribution P{x, 6) with the help of a homo- 
dyne technique. On the other hand, the Q function has 
been measured recently in various experimental schemes, 
using well-known techniques of photo detection, together 
with related experiments of phase measurement |5| . 

Apart from cavity QED, a single trapped ion is the 
other testing ground of the intriguing features of quan- 
tum mechanics The motion of a single trapped ion 
can be easily modified using laser light, and decoherence 



in such a system can be made nearly negligible during 
long times. Using these properties, several proposals have 
emerged dealing with the preparation of non-classical 
states of motion. Just recently, the first observation of 
non-classical states such as Fock states, and squeezed 
states has been reported 0. Hence, the next step of re- 
search is to characterize these states. Given the analogy 
between cavity QED and a trapped ion interacting with a 
laser, one could imagine that some techniques developed 
in the framework of cavity QED can be immediately tran- 
scribed to the ion system. For example, one can char- 
acterize the motional state by measuring the evolution 
of the ion population inversion [p|-pl| . Endoscopic tech- 
niques, for example, permit a complete state detection if 
there is no statistical uncertainty in the state prepara- 
tion process. Unfortunately, this method does not allow 
to recover the whole density matrix describing the ion 
motion. Moreover, the mentioned analogy is only valid 
in the Lamb-Dicke regime, whereby the motion of the 
ion is restricted to a region smaller than a wavelength, 
which limits the applicability of these methods. Thus, it 
would be desirable to have a method to recover the full 
information about the motional state of an ion valid for 
more general situations. 

In the present letter we propose a novel realization of a 
phase space tomography to determine the motional state 
of an ion in a harmonic trapping potential. In contrast to 
a recent proposal ||l2| that addressed the same question, 
our scheme is not restricted to the Lamb-Dicke regime 
and can be extended easily to more than one spatial di- 
mension. Furthermore, an implementation of this idea 
is feasible with present experimental setups. Specifically, 
we will present procedures: 

(i) to measure the Q{a) function, and 

(ii) the quadrature distribution function P{x,9). 

Our model consists of a single ion, trapped in a har- 
monic potential oscillating with a frequency v. The in- 
ternal structure of the ion will be specified later in the 
context of the measurement of the motional state. We 
use a density operator p to describe this unknown state 
of the particle and represent it in the Fock basis of the 
harmonic oscillator, i.e. 

oo 

P= ^ PnryiW) {m\. (2) 
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Let us first show how the Q function given in Eq. ( |la| ) 
can be measured experimentally. For this purpose, we 
re-express the Q function as 



Q{a) 



■(0|p|0), 



where 



p = u{\aie)pu\\aie). 



(3) 
(4) 



Here, C/^(|a|,6') = 7^t(0)P(|a|) is the unitary trans- 
formation that is given by the displacement operator 
'D{a) = exp(aa''' — a*a) and the phase shifting operator 
H-iO) — exp(— i^a^a) that acting on the vacuum create a 
coherent state | jal e*^) = C/^(|a|, 6*) |0). As usual, and 
a denote creation and annihilation operators that obey 



[a, a 



1. 



According to Eq. (||), one has to determine the prob- 
ability of the state represented by p to be in the ground 
state of the harmonic potential, in order to measure this 
Q function. Note that such a state is related to the orig- 
inal p by the unitary transformation C/(|ck|,0). Conse- 
quently, the identification of this transformation with a 
physical process would enable us to measure the Q func- 
tion. In the context of an ion trapped in a harmonic 
potential, this identification is as follows. The operator 
7?. (6') corresponds to the free evolution, whereas the op- 
erator I?(a) corresponds to a sudden displacement of the 
harmonic trap. A pictorial representation of these opera- 
tions in phase space is shown in Fig. 1(a). Thus, in order 
to measure the Q function in a trap, one simply has to 
perform the following steps: (i) Wait a particular time t 
while the ion evolves freely in the trap. This gives it the 
appropriate phase shift according to 9 = vt. (ii) Sud- 
denly displace the center of the trap to the right for a 
distance d, so that |q| = \f^d. (iii) Finally, measure 
the probability of the ion to be in the lowest motional 
state |0). 

To achieve this last step, one may use the internal 
structure of the ion. Typically, it consists of three levels 
\g), |e), and \r), where \g) — > \r) is a dipole-forbidden 
transition or Raman transition, whereas \g) — > |e) is a 
dipole-allowed transition. Initially, the ion is in the in- 
ternal ground state \g). After the step (ii) a laser beam 
is tuned to the lower side-band of the \g) \r) internal 
transition. One can then transfer completely the popu- 
lation of the ground states \n,g) (with n = 1,2, . . .) to 
the excited states \n — l,r) coherently as described in 
Rcf. ||l5| by an adiabatic sweep of the laser frequency. 
After this population transfer one can switch on another 
laser, this time on resonance with the transition \g) — > |e) 
[see Fig. 1(c)]. The appeareance of fluorescence indicates 
the presence of population in the \0,g) state. One can 
repeat the same sequence of steps in order to determine 
the probability of the ion to be in the ground state. An 
alternative (more sophisticated) way of measuring this 
probability may be achieved by detecting the collapses 
and revivals in the population inversion, since this tech- 
nique provides the whole population of the Fock states 



Up to now, we have not addressed the question of the 
final reconstruction of our state from the experimental 
data of the Q function that is obtained in this manner. 
In principle, one could use the method ||l^ that relates 
this function to the matrix elements Pn,m- This is, how- 
ever, impractical, since it requires the n-th and m-th 
derivatives of the Q function, i.e. their knowledge over 
a continuous interval. Another possibility would be to 
assume that pn,m = for n,m > rimax, for a given nmax- 



The measurement of Q{ai) for 



(independent) val- 



ues of Qfi would allow us to find pn,m by simple matrix 
inversion. This procedure is also of limited usefulness 
since small deviations from the exact values of the Q 
function (such as experimental uncertainties) give large 
errors in the reconstruction. This is due to the fact that 
Q is the smoothest function of all s-parameterized quasi- 
distributions. 

An alternative way of reconstructing the state of a 
quantum system is by means of quantum tomography. 
Tomography is an experimental tool used in several ar- 
eas of research which allow us to reconstruct an un- 
known object from measured data. In the context of 
quantum physics, the data we will have to measure are 
the so-called quadrature distribution functions given by 
Eq. (lb) where \x,9) = TZ^6)\x) is the eigenstate of the 
operator x{9) — TZ^ {9)xTl{9), with eigenvalue x {x is the 
dimensionless position operator of the harmonic oscilla- 
tor mv/h 1). This distribution is equivalent to that 
given by the marginal distribution for quadrature compo- 
nents using the Wigner function description of the state 
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Our scheme for the measurement of P{x,9) is based 
on the well-known property of the squeezed states 



\x,9) = lim J\ft\a, 

|e|^oo 



(5) 



exp[(e*a2 _ eat )/2] 



and a = xe'^/V^ 



where |a,e) = V{a)S{e)\0), 5(e 
is the "squeeze" operator, e = |e|e 
(0 < 6* < tt). As proper position eigenstates are 
not normalizeable, there is a constant of proportionality 
J\fe — (exp (2|e|)/(47r))4 that increases with the degree of 
squeezing. As before, we can use these states to reexpress 
the quadrature distribution in the form 



P{x,9) = lim |AA,p(0|p|0). 



(6) 



Here 



p^U{\e\,\a\,9)pU\\el\al9), (7) 



where Wi\e\,\a\,9) = n^9)V{\a\)S{\e\) denotes the op- 
eration that creates a squeezed state and we used fur- 
thermore the property S{e)= n^9)S{\e\)n{9). Thus, to 
measure the quadrature distribution one has to find the 
physical processes that correspond to the unitary oper- 
ators TZ, V and S. The first two are the same as those 
needed for the experimental determination of the Q func- 
tion. On the other hand, it is well known that sudden 
changes in the frequency of a harmonic oscillator lead to 
squeezed states, a process that can be readily achieved 
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in a trap, just opening or closing the harmonic potential 
1^^ . In particular, changing the trap frequency from v 
to v' leads to a squeezing parameter |e| = ^ln^. 

Thus, in order to measure the quadrature distribution 
in our trap one has to follow these four steps: (i) Wait 
for a time i, such as 6* = vt. (ii) Perform a sudden dis- 
placement of the center of the trap to the right a distance 
d, so that |a| = yj^d. (iii) Change the trap frequency 
instantaneously from v Xo v' . (iv) Determine the popu- 
lation of the motional ground state. Note that the steps 
(i,ii) and (iv) are the same as before. 

We are now in position to extract the full informa- 
tion about the unknown quantum state starting from the 
quadrature distribution. This can be done as in the case 
where one measures the quantum state of light by means 
of balanced homodyne detection. As has been shown 
by Vogel an Risken one can reconstruct the Wigner 
function W^(a;, p) by means of the so-called inverse Radon 
transformation. Alternatively, one can use one of the al- 
gorithms that have been developed for reconstructing the 
density matrix directly from discrete measured data p9[ . 

To illustrate this procedure, we have simulated numeri- 
cally the reconstruction of a quantum state. We assumed 
the system is prepared initially in a " Schrodinger-cat" 
state of the form 

where |a) is a coherent state. This is a highly non- 
classical state, and can be easily produced in the trapped 
ion system |Q . In Fig. 2 we plot the real part of the re- 
constructed matrix elements Pn.m corresponding to the 
initial state Eq. (^) with a — 1.5. We have taken the 
values of P{x,9) for a set of points {xi,9j) and with a 
finite e. Starting from these data we have reconstructed 
the state pn,m using the algorithm of Leonhardt et al 
po| . We have selected a uniform grid of points corre- 
sponding to values of x ranging between ±4[?i/(mj/)]^/^, 
and a uniform grid of Ng points for < < tt. Figure 
2(b) corresponds to a squeezing parameter |e| = 2 and a 
grid Nx X Ng = 30 x 30 points. The reconstructed state is 
indistinguishable from the original one. We have checked 
that even for |e| = 1 the obtained state is remarkably sim- 
ilar to the original one. We have also tested the depen- 
dence of the reconstruction on the number of grid points. 
In Fig. 2(b) we have taken a grid of TV^, x TV^ = 30 x 15 
points, keeping the squeezing parameter |e| = 2. In this 
case, the reconstruction is also quite faithful. Reducing 
the number of grid points causes small residual back- 
ground structures. On the other hand, reconstructing 
density operators that involve higher Fock states (in- 
creasing a) requires an increased range of x values and 
a larger number of grid points, since it is necessary to 
resolve the oscillatory behavior of these states. Finally, 
in a real experiment one cannot measure the probabil- 
ity distribution P{xi, 9j) with arbitrary precision, due to 
the fact that the number of measurements is always fi- 
nite. We have simulated the statistical error caused by 



the finite sampling number by truncating the values of 
P{xi,9j) to one decimal digit. That is, we have approx- 
imated each of the exact values by one of the following 
numbers 0.0, 0.1, . . . , 1.0. The results of this simulation 
are shown in Fig. 2(c). In this case, the grid size is again 
Nx X Ng — 30 X 30 points, and |e| = 2. The reconstruc- 
tion still resembles the original one even in the presence 
of these uncertainties. Therefore, it would be enough to 
perform about 100 measurements per grid point to ob- 
tain the density matrix. Obviously, states with a larger 
phonon number will require more measurements. 

In summary, we have presented a scheme to measure 
the quantum state of motion for a single ion confined 
in a harmonic potential (statistical mixtures as well as 
pure states). It is based on the detection of the ground 
state population of the trap after a sudden change of 
the trapping potential. We wish to emphasize that the 
effect of a sudden displacement of the trap center and 
a sudden opening of the trap can be obtained (in prin- 
ciple) by a single non-instantaneous process that yields 
the same symplectic phase space transformation. More- 
over, these two operations can be mimicked using Raman 
pulses with two lasers of frequencies differing by p and 
2 J/, respectively |0,|l^. Note, however, that in this case 
the scheme only works in the Lamb-Dicke limit. Finally, 
it is straightforward to generalize the schemes presented 
here to measure the quantum state of motion in two and 
three spatial dimensions. This can be done by moving 
and opening the trapping potential along different direc- 
tions. 
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FIG. 1. Phase space representation of the operations 
(phase shifting, displacement and squeezing) required to mea- 
sure the Q{a) function (a) and the quadrature distribu- 
tion (b); (c) Level scheme and laser configuration for the de- 
tection of the trap ground state population. 

FIG. 2. Real part of the reconstructed density matrix 
elements Pn^^ for a "Schrodinger-cat" state with a = 1.5; 

(a) X Ne = 30 x 30 and |e| = 2 ; (b) iV;^ x iVe = 30 x 15 
and jej = 2; (c) same parameters as in (a), but with P{x,d) 
rounded to one decimal digit. 
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Fig.2a 




Fig.2b 




Fig.2c 




